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Abstract 

We discuss stochastic derivations, stochastic Hamihonians and the 
flows that they generate, algebraic fluctuaion-dissipation theorems, etc., 
in a language common to both classical and quantum algebras. It is con- 
venient to define distinct notions of time-ordered exponentials to take 
account of the breakdown of the Leibniz rule in the Ito calculus. We 
introduce a notion of quantum Stratonovich calculus and show how it re- 
lates to Stratonovich-Dyson time ordered exponentials. We then use it to 
demonstrate a natural way to add stochastic derivations. 



1 Introduction 

Symmetries play a considerable role in Mathematical Physics, particularly when 
determining dynamical flows with required invariant properties. Much of our 
insight comes from being able to go from the infinitesimal generators (having 
desirable symmetry features) to the flows themselves. This intuition requires, 
to a large extent, the generators to be derivations. When we consider stochastic 
flows, the ltd calculus looses the Leibniz rule and with it much of our physical 
intuition. This is why the Stratonovich calculus PP is usually preferred by physi- 
cists. In addition, the Wong-Zakai approximation theory tells us that we can 
approximate a stochastic flow with a differentiable random flow: if the approx- 
imations have some symmetry property then the limit Stratonovich equations 
will retain this property before it is typically lost in the conversion to Ito form. 

The aim of this article is describe algebraic notions such as derivations, 
their stochastic analogues, dissipation, etc., in a common language that applies 
to classical and quantum systems. Whereas there have been many attempts to 
construct general classical analogues to quantum stochastic flows, the spirit of 
the quantum flows is arguably best captured by stochastic flows on symplectic 
manifolds preserving the Poisson brackets. This was investigated by Sinha 
and we develop somewhat the algebraic similarities. Given the suitability of C*- 
algebras for modelling quantum mechanical variables it is natural to consider 
algebras of functions on Poisson manifolds as an intermediary between the com- 
mutative and noncommutative cases. This view is strengthened considerably by 
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the deep analogies known to exist between Poisson and operator algebras [3], 
It is natural to return to these analogies when considering stochastic flows 
describing irreversible dynamical evolutions. 

One of the most notable omissions from quantum probability has been a 
general theory of a quantum Stratonovich calculus. This is largely due to the 
bias, already prevalent in classical probability, towards the Ito calculus on the 
grounds that of martingale stability. Nevertheless, in problems, such as stochas- 
tic process on manifolds [S] [S], the Stratonovich calculus is used in on entirely 
equivalent footing as the Ito calculus and, because it does maintain the Leibniz 
rule, allows one to see the underlying differential geometric structure. Mo- 
tivated to develop these idea for algebras of quantum observables, the author 
considered Stratonovich integral for quantum diffusions [7|, however, in the later 
development of this ideas, much of the physical insight came from limit theo- 
rems - weak coupling and low density - and this lead to a formulation in terms 
of quantum white noises: see [H], [HI- Formally, the familiar features, like a 
Hamiltonian nature for the flow, are evident in the Weyl ordered form (the ana- 
logue of Stratonovich version) while the Wick ordering distorts these features 
and produced the analogue of the Ito version. While quantum white noises of- 
fer the most intuitive approach, and should amenable to a rigorous treatment 
in some extension of the Hida theory, they have not found favour with either 
the quantum probability or the mathematical physics community. However, it 
is possible to give a sufficient account of events from within the quantum Ito 
calculus, and this is what we address in this paper. 

With regards to applications, it is convenient to distinguish three different 
notions of time-ordered exponentials when exponentiating quantum stochastic 
integrals. The usual Dyson form, with Ito differentials, does not exponentiate a 
derivation-valued process to a homomorphic map. Instead we must either use 
Stratonovich differentials or exponentiate over time steps. The former strategy 
was applied using quantum white noises in a series of papers by the author [S], 
[TU). while the latter was investigated in jJT] and These alternatives to the 
Ito-Dyson form coincide in the special case of quantum diffusions. 

Let 21 be the *-algebra of operators modelling a system and its environment. 
The space of all linear maps L :by 21 1-^ 21, having the reality property L {X'^^ = 
L {X)\ will be denoted by C (2t). The dissipation of such a map is defined to 
be the bilinear mapping : 21 x 21 21 given by 

'Z)l{XX) = L{XY)-{LX)Y - X{LY). (1) 

li Dl is zero then L is called a derivation, li v E C (21) is a derivation on the 
algebra if we have the Leibniz property 

V {XY) ^v{X)Y + Xv {Y) , 

for all X, y €E 21. Let us introduce the notations 

on tv / \ \ ^ on / \ 

n times n— 
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then {XY) = J2m (m)^°" (^) (^) and 

e*" (Xy) = e*" (X) Y + Xe*" (Y) . (2) 

Therefore, derivations act as the generators of flow maps that preserve the 
algebraic structure (homomorphisms) . 

We introduce the following causal structure: for each f > 0, let Stj] be the 
*-subalgebra of 21 modelling the system and its environment up to time t, then 
we assume that we have the isotony condition 21,,] C 2lt] whenever s < t. The 
family {Slj] : t > O} is called a filtration. Let E^j : 21 21^] be a projective 
conditional expectation 

A flow is a family {^t,s ■ t > s > 0} of maps in £ (21) with the properties 

i) ^t^r ° ^r,s = ^t,s, whenever t > r > s; 

ii) limtxs ^t.s = id, the identity map on 21. 

The flow is said to be adapted to the filtration if ^t,s {^s]) ^t] whenever 
s < t and is said to be homomorphic if $t,s i^Y) = (X) (Y) for all 
X,Y€^s]- 

Given a flow, we define the maps {i(,s : t > s > 0} of maps in C (21) by 
Lf^s = ^t,s — id. Let us write dL^ for Lt+cit,t and we may think of this as an 
C (2t)-valued measure which we refer to as a stochastic generator. Given two 
such measures dL^ and dL'^, we define their mutual quadratic variation over 
time interval [S, T] as 

/ dLt o dLi = lim >^ L (t-j+i +ti,ti) o L' (f,+i + tjAi) 

where we have the limit over all partitions P {S, T) of [S, T] into sub-intervals of 
maximum length \P {S, T)\. In general, for a family {-K't,s : t > s > 0} of maps 
in £ (21), if we have that, for every finite interval [S, T], 

lim K (ti+i + ti,ti) = 

then we write dKr = o{dt). We say that a flow is regular, or deterministic, 
if dLt o dLt = o (dt). In general, however, this docs not hold and we typically 
encounter the Ito rule for differentials of compositions, viz. 

d{Lo L') = (dL) oL' + Lo {dL') + (dL) o {dL') 

with the last term, the Ito correction, being non-zero. (Here "=" means equal 
up to terms of order o (dt).) 

A flow $ is said to admit a forward derivative (•) e £ (21) if the following 
limits exist 

ut (X) := lim ifit] [^t+r,tX - X] (3) 

T— >0+ T 
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for each t > 0, X S 2lt] . The forward derivative will generally not be a derivation! 
If the flow admits a forward derivative ut then a difference martingale M.t.to (■) £ 
C (Sljg]) is defined by 

iX) = X+ f {X) ds + Mt.to {X) , {to < t) . 
For s < r < t, Mt,s = Mt,r + Mr,s- 
1.1 Ito-Dyson Exponentials 

The flow can be reconstructed from dLt. We have the (Ito) differential equations 

d^t^to (•) = dLt o <i>t,to (•) (4) 
with initial condition Ihntito ^t.to = id. We write 

$t,to =f/i3exp|^ dlj (5) 

where Tjd is the Dyson time-ordering symbol and we refer to right hand sides 
as Ito-Dyson time- ordered exponentials. They can be developed as an, at least 
formal, Picard series 

$t,t„ {X) = X+ f dLt, {X) + f dLt, ( f ' dLt, {X) 

Jto Jto \Jto 

oo ^ 

= '^^t- ° • • • ° '^^'i ■ (6) 

Here A„ (a, b) denotes the simplex consisting of all n— tuples {tn, - ■ ■ ,ti) ordered 
so that b > tn > ■ ■ ■ > tl > a. So we have, again formally, 

TiD exp M dLV dLt„ o ■ • ■ o dLt, . 

I "'to J „^o"'^"(*0'*) 

Because the Ito rule means a breakdown of the Leibniz identity, it typically 
means that if dLt {■) is equivalent to a derivation, but is not regular, then the 
flow $ we construct will fail to be homomorphic. Indeed, the condition on dLt in 
order for the flow to be a family homomorphisms is that its dissipation balances 
its fluctuations ^31) '^i^- 

SdL. = {dLt) ® {dLt) ■ (7) 
More explicitly, dLt {XY) ^ dLt{X)Y + XdLt{Y) + dLt{X)dLt{Y) . 
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1.2 Exponentiated Dyson Exponentials 

We may refer to ^t.to = e*^*"*")" as the autonomous flow generated by velocity 
field V. More generally, we could take a family {vt : t > 0} of derivations and 
consider the flow generated by dLf (•) = Vt (•) dt, that is, Vt (•) is the instan- 
taneous velocity field. This time, we can write the non-autonomous flow as 
^t,to = T/D exp I J^^ Vgds^. We observe that 

^t+dt,t (•) = id (•) + Vt (•) dt = e<-^'^' + o {dt) . 

Unfortunately this cannot be the case when we consider differential genera- 
tors dLt that are not regular. The reason is that, because of the Ito calculus, 
compositions {dL^)"" need not be o{dt) for n>2. Let us write in general 

dHt = e"^^* -id=Y, {dLtY"" (8) 

and define the exponentiated Dyson, or i^D-type, time-ordered exponential to 
be Ted exp | J^* dL^ = Tid exp | J^* rfijj .This time, if the dLt are derivations 
then the e**^* (•) will behave as homomorphisms, and so too will 

Ted { | = T/d exp | ^ (e'^^ - i d) | . (9) 

1.3 Stratonovich-Dyson Exponentials 

As we have seen above, the Ito calculus implies that the Dyson time-ordered 

exponential T/£)exp | J^*^ rfij will not generally yield a homomorphism if the 

dLt (•) are derivations. There is another strategy for producing homomorphisms 
and that is to use the Stratonovich calculus instead of the Ito one. The Ito rule 
(for ordinary products) is that d(XtFt) = Xt {dYt) + (dXt) Yt + {dXt) (dYt). The 
Leibniz rule may be formally restored as d {XiYi) = * (dYj) + {dXi)*Yt where 
we introduce Stratonovich differentials Xt * {dYt) = Xt {dYt) + \ {dXt) {dYt). 
For compositions of infinitesimal generators we have the similar relation 

d {X, o Y,) = Xt o {dYt) + {dXt) o Yt + {dXt) o {dYt) , 

with the Leibniz rule salvaged as 

d {Xt o Yt) =Xt® {dYt) + {dXt) ® Yt 

where now Xt ® {dYt) =Xto {dYt) + \ {dXt) o {dYt) and {dXt) ® = {dXt) o 
Yt^\{dXt)o{dYt). 

We define the Stratonovich-Dyson, or SD-type, time-ordered exponential 

^t,to = TsD exp I dLg I to be the solution to the integro-differential equation 



<i>t,to {X)=X+ I dGs {^s,to {X)) 

J to 
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where dGt is the Stratonovich differential 

dGt o $ = dLt (g $ = dLt o $ + ^dLt o dLt o $. 

that is 

rfGt =dLt + ^ {dLt)°^ . (10) 

We have now defined $t,to = T^SD exp | J^*^ rfLgj as the solution to 

^t+dtM (■) = d^t ® ^t,to (•) = dGt ° *t,to (•) 
with \imtito ^t,to = '''d- As the Leibniz rule is observed in the Stratonovich calcu- 
lus, TsD exp I J^* (iLsj will be a homomorphism when the dLt are derivations. 

1.4 Properties 

If the flow is regular, say dLi (•) = Vc (•) (if for some Vc G /C (21) and {dLt)°"' = 
for n > 2, then the various notion of timc-ordcrcd exponentials coincide. We 
say that the flow is a diffusion if (dL^)"'^ is equivalent to At (■) dt for some 
At ^ C (21), however, {dLi)°^^ = for n > 3. In this case, we have the truncation 

dHt = e''^* -id = dLt + ^ [dLtf^ 
and so dHt = dGt and therefore 

Ted exp < / dLg I = Tsd exp < / dL, 

This is however a fluke which cannot be expected to hold for stochastic processes 
other than diffusions. 

The ED-type exponential occurs when we consider discrete time-approximations. 
We have for instance 

Ted expi dL\ = lim exp {^tN,tN-i } o • • • o exp {$ti,to} 

Uto ) max|t3 + i-tj|->0 

= lim exp < / dLg > o • • • o exp < / dLg > 

maxltj + i-tjI^O J IJto J 

where the limit is over all partitions t = tN > In-i > ■ ■ ■ > h > to- Let us 
suppose that the stochastic generator has essentially commutative increments, 
that is, 

whenever the intervals (s, t) and (s', t') do not overlap, then 
TBnexpj^ dL^j = exp^t,to- 
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The iSD-type exponential occurs when we wish to approximate stochastic 
flows by regular ones. Let v^:^^ (•) be velocity fields parametrized by i > and 
A > 0. Suppose that we have the limit 

hm f\i^H-)ds= fdLA-) 

^^UJto Jto 

where dLt is possibly stochastic. We may then expect the limit 



lini f exp 1^ vi^Us^ = TsD exp |^ di^j ■ 

Evidently we should not have expected the Dyson exponential on the right 
hand side as limit should be a homomorphism. In the theory of approximating 
stochastic differential equations by ordinary differential equations, it is well 
known that the Stratonovich calculus is the one that best anticipates the limit 
form. 



1.5 Examples from Mathematical Physics. 

Classical mechanics Let 2t be the C°° functions on a manifold M. Here 
21 is a commutative algebra (though not generally a C*-algebra) with respect 
to pointwise multiplication, and derivations correspond to the tangent vector 
fields. In this case, the dissipation is known by several different names: the 
Gamma operator, I'operateur carre du champ, the cometric operator, etc. (see 
Meyer's appendix to [H]). 

If, however, M is also endowed with a Poisson brackets which is, of course, 
a bilinear mapping. We may take our product to be the anti-symmetric, non- 
associative one given by the choice f-kg = {/, g}. In this case the 7k-automorphisms 
are those maps that preserve the Poisson brackets. Likewise S^L will now de- 
termine the extent to which a semi-group e*^ destroys the Poisson structure: 
explicitly, we have Dl {f,g) = L{{f,g}) ~ {Lf,g} - {f,Lg}. Given a real 
function /i G 2t, the Hamiltonian vector field Xh generated by h is defined by 

if) '■= if, h}. By the Jacobi property of Poisson brackets, Xh will be a 
Poisson-derivation 

Xh{f,g}^{Xhf,g} + {f,Xhg}; 

and by the Leibniz property of Poisson brackets it will also be a tangent vector 
field. The Poisson manifold is said to be Poisson-simple if the only maps having 
Poisson-dissipation zero are the Hamiltonian vector fields For symplectic 
manifolds, all Poisson derivations are locally Hamiltonian generated. 



Quantum Mechanics Next let 21 be an algebra of operators acting on a 
Hilbert space Ti., with the taking of adjoints f as the usual involution. The 
operator product is non-commutative and to a certain extent carries out the 
role played by both pointwise multiplication and Poisson brackets in classical 
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mechanics. The hnear maps on 21 i— > 21 are sometimes referred to as super- 
operators. For a self-adjoint operator /i £ 21, we define the super-operator Xhby 
Xfi if ) J [/, h]. If 21 is a W*-algebra, or a unital simple C*- algebra, then 
it is well-known that all real derivations take this Hamiltonian form Jl]. The 
dissipation of a super-operator, with respect to the operator product, was intro- 
duced by Lindblad |15| as a key ingredient in analyzing generators of completely 
positive semi-groups. 



2 Quantum Stochastic Flows 

Let S) = r (^L^ (R+,(it)) be the Fock space over square-integrable functions of 
positive time and let Sjt = ^ {L^ ([0, t] , dt)). The family {S)t ■ t > 0} then gives 
a filtration of Hilbert subspaces of Sj. Fixing an initial Hilbert space f), we 
consider the filtration of2l = Q3(()(g)i3) specified by 

21^] =<B([)®jOt]). 

The quartet of fundamental quantum stochastic processes may be denoted 
by ^Af^ : < > o|: these are = t (time), Al° = At (creation), Al° = A\ 
(annihilation) and A\^ = At (conservation). Their Ito table is then 

dAfdA\^ = dAf, 

with all other second order differentials vanishing. This may be written explicitly 
as HSl 



X 


dA'f 


dA 


dA 


dt 


dA^ 














dA 


dA^ 


dA 








dA 


dt 


dA 








dt 















We have the adjoint relations yA'^'^j = ^f": that is. At is self-adjoint and A] 

is indeed the adjoint of At We remark that the combinations Qt — At + A] and 
A^f = At -I- ^t + ^t + i give representations for the Wiener and Poisson processes 
respectively when we specify the Fock vacuum as state. 

A closed evolution restricted to 2lo] — S (()) may be described by the family 
of unitaries ^ 

Ut,to = T/£) exp |-*^ Hsds^ 

where {Hs : s > 0} is a family of self-adjoint operators forming what we usually 
call a time-dependent Hamiltonian. 
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Our aim is to construct unitary quantum stochastic processes on 21 of the 
form ^ ^ 

C/t.to =T/i3exp|-i^ dG,| =Tsz3Cxp|-^^ dEs^ (11) 

where 

dGt = Gap ® dAf'^ = Goo + Gio + Goi ® dAt + Gn dA* 
dEt = Eaf3 ^ dA"'^ = Eoo dt + Ew ® dAl + Eoi dAt + Ell ® dAf 

(We use a summation convention that repeated Greek indices are summed over 
values and 1.) We take the coefficients Gap and Eap to be bounded operators 
on f). Special choices of the coefficients should lead to stochastic evolutions 
driving by either Wiener or Poisson Noise, however, it is known that classical 
stochastic processes do not account for all the quantum stochastic evolutions 
we would wish to consider, and so we work with all four fundamental processes. 



2.1 Quantum Stratonovich Calculus 

We do not aim a full generalization of the Stratonovich prescription to quantum 
stochastic calculus, however, we shall give the algebraic rules which tell us how to 
transform certain Ito integrals into Stratonovich ones in a manner that parallels 
the classical theory. There are two main surprises; the first being that it can 
be done; the second being that there is an ambiguity in the definition. 

Let {Xt : i > 0} be a family of operators (a quantum stochastic process) on 
i)iS)Sj. If we have that Xt acts nontrivially only on the subspacef) Sjt] then we 
say that the process is adapted. Our definition of Stratonovich differentials will 
amount to 

{dXt)*Xt - {dXt)Yt + K{dXt)idYt), 

Xt*{dYt) = XtidYt) + n*idXt){dYt), (12) 

where k is a complex number with Ken — — . Evidently we have (dXt) * Xt + 

Xt * (dYt) = d (XtYt) from the quantum Ito formula ^21- The fact that we may 
choose the imaginary part of k means that we have a "gauge freedom" which in 
many respects is similar to that in the Tomita-Takesaki theory. The symmetric 

choice would be k = — , however, its explanation is as a damping constant, and 

physically this may be complex. 

The stochastic Schrodinger equation from above is then (ignoring the to 
dependence) 

dUt = -I {dGt) Ut = -I {dEt) * Ut (13) 

and we would like to determine how to transform between Ito and Stratonovich 
forms. Let us begin by computing (dAt^j *Ut ^ (dA"^j Ut + n (dA'^^j {dUt). 
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We have 

(dAf){dUt) = -i(dAf){G^,^dAnUt 

= -i6i0{Gi,®dA^'')Ut 

from which we see that 

(dAt) *Ut = {l- ikGh) (dAt) Ut - inG^ {dt) Ut 
_ {dt)*Ut = {dt)Ut 

and inversely 



f (dt) Ut = {dt) * U, 
(dAt) Ut = (dA\^ * Ut 
{dAt)Ut = {l-iKGii)-^ 
{dAt)Ut = il-iKGn)~^ 



(dAt) * Ut + mGio (dt) Ut] 
{dAt) * Ut + «Gio Ut 



We can next of all read off the relationship between the Ga^ and the E„/3 by 
comparing coefficients: 



or inversely 



Ell 
Eio 
Eoi 



Gn 

Gio 
Goi 
Goo 



1 — iKGll 

1 — iAcG^i 
Goi 



Gio 



1 — zkGh 



Eoo = Goo + i«;Goi 



1 



1 — mGii 



Gio 



-11 



1 + iVcEii 



-10 



-01 



1 + zkEii 



= Eoo + ii'i^i 



1 



01: 



:iO 



1 + mEii 

The relationship between the Ito and Stratonovich coefficients may be written 
more compactly using the following remarkable formulas 



Ea/3 



Ga/3 + i«;Gal: 



Gl/3, 



1 — mGii 

Ga/3 = - mEai— — — — El/3 

l-|-^«;Ell 



(14) 
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showing a duality between them. In particular the operators 1 — inGn and 
1 + jkEh are required to be invertible and inverse to each other. 

The conditions for unitarity are the isometric and co-isomctric properties 



1 = UtUj. In differential terms this is 



d ( U}Ut 



{dU}) 



(dUl) Ut + Ul (dUt) + (dC/t j (dUt) 
*Ut + ul * [dUt) (15) 

with = d {UtUf^ similarly. From the Stratonovich form dUt = —i (dEt) * Ut 

we see that it is enough to ask that {dEt)'' — dEt- The unitarity of the process 
should come down to the conditions 

(E,;,)^ = Ep^. (16) 

The process Ut^to = 'TsD^xp^—i J^*^ dEg^ is evidently unitary because we are 
using the Stratonovich calculus to time-order explicitly unitary components. If 
we use the quantum Ito calculus with dUt — —i (dGt) Ut we see that 



= d{ UjUt 1 = U, 
which implies 

= 



{ulUt) = Uj It (dGl) - I {dGt) - \ {dGt)^ {dGt) 



Ut 



i(dG\) -i{dGt)-\{dGt)UdGt) 



a/3 
t ■ 



iG^^Gi^ guarantee 



It turns out that the four equations = iG^^ — iGap ^ 
both the isometric and co-isometric properties, and therefore unitarity. It is 
well-known that the general solution to this equation is [TB| 



1. 



H -i-K^K 
2 

K 



Goo 
Gio 

Goi = KtW 
Gn = z(W-l) 

with W unitary, H self-adjoint, and K bounded but otherwise arbitrary. One 
can readily check that the coefficients Gaf) will satisfy these conditions once 
(Eq,^)^ = E^q, with the explicit choices 

1 — iK*Eii 



W = 
K = 
H = 



1 + zkEii 
1 



1 + 
Eoo 
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+ Im/tEoi 



1 



1 + zkEh 
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So, although the stochastic generator dGt does not look self-adjoint, the process 
Ut^to = T/£) exp |— i Jj*^ (iGs| is nevertheless unitary. 

The flow is then given by the family of maps 

<^tM {X) = Ul, {X) Ut,t, 

which defines a system of homomorphic flows on 21. The quantum stochastic 
differential equation for ^t,to (X) is 

d$t,to = dLt o $t,t„ 

with 

dL,i-) = £ip {■)^dAf 

where C^p (X) = z (g*,„)^ X-zX (g^^) -i {GlJ X (g^^) and G^^ = (G^^). 
The stochastic derivation property is then stated in the form 

that is, Ca,3 {XY) - £„/3 {X) Y - XC^p {Y) = L^i {X) Cip {Y) . These are the 
well-known structure equations for non-commutative flows |17) . 



2.2 Remarks 

In principle, the algebraic manipulations can be extended to time-dependent 
E„/3 G 05 (f)) and more generally to adapted coefficients. We may also generalize 
to N Bose noises A"^ where now the Greek indices run over 0, 1, • ■ • TV and this 
leads to a tensorial version of the equations (14). In this case we are free to 
introduce additional gauge degrees of freedom. 

We mention that we have the following approximations theorem which is the 
quantum analogue of the Wong-Zakai result and which justifies our construction 
of quantum Stratonovich calculus. Let a\ (A) be Bose fields on a Fock space io^'' 
correspond to some physical system which we shall refer to as the reservoir, we 
have chosen to parameterize them by time t and also a scale parameter A > 0. 
For A fixed, we consider canonical commutation relations of the type 

[at{\),a\{\)\^Gx{t~s) (17) 

where G\ is the two point function and is assumed to be a regular function of 
the time difference. We shall assumed that G\ (•) is integrable with G\ — 1 
and we naturally require that 

Ga (-r) = Ga (r)* . 
Let K = Ga then Re k = ^ and Ga — k* . 
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We now consider what happens if, in the Hmit A — > 0, we have 

G,{t)-^6{t) (18) 

in the sense of Schwartz distribution. In particular, we wish to study the asymp- 
totic behavior of the unitary Ut (A) coupUng a given system with state space f) 
to the reservoir with an interaction Hamiltonian Tf (A). Here Ut (A) is given as 
the solution to ^ 

Ut{\)^l-i ( T, (A) Us (A) ds. 



Theorem 1 For the interactions on f) ® of the type 

Tf (A) = Ell ® a\ (A) at (A) + Eio ® a\ (A) + Eqi ® at (A) + Eoo ® 1 (19) 

with Eap e /8([))j Ell (ind Eqo self-adjoint, Eio = Eqi, and ||kEii|| < 1, the weak 
matrix limit of Ut (A) is described by the unitary quantum stochastic process Ut 
given by the Stratonovich-Dyson time-ordered exponential 

Ut = TsDCxp|-i / dEs 



to 



with dEt = Ea0 <E) dA"^ . 



The condition ||kEii|| < 1 is required to ensure that multiple scatterings 
diminish rather than augment amplitudes - it also allows (14) to be expanded 
in a geometric series. A simpler version of this result, applicable for commuting 
coefficients only, was given in ^5]. The convergence also applies to the Heisen- 
berg dynamics and so we get convergence of the regular pre-limit flow to the 
stochastic quantum flow. A similar set of formula arise when the a\ (A) are 
replaced by Fermion fields the fundamental processes A'^^ now being the 
Fermi analogues. 

2.3 Addition Rules for Stochastic Derivations 

Let I fiLf"'' [ be a finite collection of stochastic derivations: their sum is not 

typically another stochastic derivation. In general, ^„ dL^^^ + dF defines a 
stochastic derivation only for some suitable choice of "Ito" correction F. For 
the quantum problem, we realize each stochastic derivation dLt as a function of 
the operators E^/j, i.e. dLt = dGt {^ap)- 

The natural procedure is then to consider the total stochastic Hamiltonian 
dEt = rf-Et"^ = X^n ^^ap ® dA^^ . The corresponding stochastic derivation 
is then 
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2.3.1 Examples 

i) When the E^"-* = 0, the relation is simply 

K = ^K("); H = ^H(") -ReK^K(")tK(") 



ii) (For simplicity, take k = i.) Let W'^^^W'^^ be commutative unitaries 

W<°'-1 

w(°)+i ■ 



related to E^"-* and E^^j^^ by the preceding relations That is, E^"-* — 2i^jj^^y^. The 



l-z(E<t'+E<';')/2 

composite miitary is W = — — (TH — which, after some algebra, becomes 

l+i( Ejj +Ej^j^ j /2 

. ^ (3 + W(") +WW - W('')W(''))^ ... 
W = W^'')-^ — ■ ■ ^-W^*") 



In 1 a formula K := L + M + [[L,M]] for the sum of two stochastic 
derivations, L and M, is given. There the bracket is [[L, AI]] is their mutual 
quadratic variation defined by 

[[L,M]] {t,dt) = dLtodMt 

where o denotes composition in C (21) . We have developed a generalization to 
arbitrary many summands. 
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